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Table 1. Revised orbital parameters of the reference model 



Parameter 



Value 



Edge radius 


142 pc 


Minimum orbital velocity 


179 km s" 1 


Maximum orbital velocity 


254 km s" 1 


Current orbital velocity 


184 km s" 1 


Apocentre 


7.65 kpc 


Pericentre 


6.10 kpc 


Current Galactocentric distance 


7.60 kpc 


Orbital eccentricity 


0.11 


Orbital phase 


+0.97 



The paper "The Tidal Tails of 47 Tucanae" was pub- 
lished in Mon. Not. R. Astron. Soc. 423, 2845 (2012). An 
error was discovered, which affected the results reported in 
that paper. This error is corrected here. 

We used, for our pr oper motion of 47 Tucana e (47 Tuc), 
the values published bv lAnderson fc K ing (2004), however, 
we incorrectly adjusted the published values of \i a by cos 8. 
These values were already published with the cos S included 
so this was unnecessary. Our new reference model for this 
erratum uses proper motions: 



[la cos S — 5.64 mas yr 

fj,s — — 2.05masyr _1 



(1) 
(2) 



All other input parameters were correct and are left un- 
changed. All references for, and reasons behind, our choices 
for the input parameters can be found in the original paper. 
The newly computed models are presented in the following. 



1 THE ORBIT 

The last 500 Myr of the (x,y) and (x,z) orbital paths of 47 
Tuc derived from the new proper motion and radial veloc- 
ity of the cluster are shown in Figure [T] Due to the signif- 
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icant change in proper motion, various orbital parameters 
have been altered. The new quantities, defined as per the 
original paper, are listed in Table [1] We particularly note 
that (for our choice of Galaxy model) the orbit of 47 Tuc 
is less eccentric than originally thought (e = 0.11) and that 
47 Tuc is very close to reaching the apocentre of its orbit 

(Port = +0.97). 



2 THE TIDAL TAILS 



Figures[2]|8]in this erratum are based on our new models, and 
correspond to Figures 3-9 from the original paper. Table [2] 
summarises our corrected results. Qualitatively, the results 
are effectively equivalent, but the shape of the tidal tails and, 
therefore, the positions of their overdensities have changed. 
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Table 2. The density and Galactic coordinates of the first and second order overdensities, the proper motion, final (present day) cluster 
mass, Heliocentric distance and escape conditions for each of our corrected models. Only the trailing tail is represented here for brevity. 
Information is for the peak density in each overdensity only. Note that the approximate average density of the stellar background in the 
regions of the first and second tidal overdensities is 3000 stars per square degree, based on 2MASS point sources. 



Model 


Si 


s 2 


(l,b)i 


(l,b) 2 


fj, a cos S 


US 


xlO 6 


D 


Esc 




* / deg 2 


*/dcg 2 


deg 


deg 


mas yr — 1 


mas yr 


M Q 


kpc 




1 


85 


115 


(-60.4,-34.5) 


(-56.7,-26.5) 


5.64 


-2.05 


1.0 


4.02 


warm 


2 


55 


70 


(-60.5,-34.5) 


(-57.0,-29.2) 


5.64 


-2.05 


1.0 


4.02 


hot 


3 


120 


110 


(-56.5,-34.9) 


(-52.2,-27.7) 


7.05 


-2.05 


1.0 


4.02 


warm 


4 


100 


100 


(-62.0,-38.7) 


(-61.6,-31.0) 


4.23 


-2.05 


1.0 


4.02 


warm 


5 


85 


110 


(-56.5,-28.0) 


(-69.2,-33.6) 


5.64 


-2.56 


1.0 


4.02 


warm 


6 


95 


130 


(-60.3,-34.0) 


(-56.9,-24.7) 


5.64 


-1.54 


1.0 


4.02 


warm 


7 


85 


120 


(-60.4,-34.3) 


(-56.5,-25.5) 


5.64 


-2.05 


1.1 


4.02 


warm 


8 


95 


130 


(-60.2,-35.0) 


(-57.0,-27.0) 


5.64 


-2.05 


0.9 


4.02 


warm 


9 


65 


85 


(-63.7,-32.6) 


(-60.4,-23.5) 


5.64 


-2.05 


1.0 


3.30 


warm 


10 


120 


125 


(-57.7,-35.7) 


(-53.7,-29.0) 


5.64 


-2.05 


1.0 


4.70 


warm 



This paper has been typeset from a TffjX/ JOTgX file prepared 
by the author. 
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Figure 1. The orbital path of 47 Tuc from 500 Myr ago to the 
present time in the (x,y) plane (top panel) and in the (x,z) plane 
(bottom panel), based on parameters from the corrected proper 
motions (see text). The present day position of 47 Tuc is marked 
as a plus (+) symbol and the position of the Sun is marked as a 
cross (x). Corresponds to Figure 2 of the original paper. 
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Figure 2. The final timestep for the reference model. Corre- 
sponds to Figure 3 from the original paper. 
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Figure 3. Same as Figure[2]except with warm (top) and hot (bot- 
tom) escape conditions. Corresponds to Figure 4 of the original 
paper. 
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Figure 4. Contour representation of Figure [3] with "warm" es- 
cape conditions in the top panel and "hot" escape conditions in 
the bottom panel. Corresponds to Figure 5 of the original paper. 



-20 
-25 
-30 
-35 
-40 
-45 
-50 
-55 

-20 
-25 
-30 
-35 
-40 
-45 
-50 
-55 



47 Tuc today 
Orbit 

|x a cos(8) = 4.23 mas/yr 
King radius 




-10 



-20 



-30 -40 
I [deg] 



-50 



-60 



-70 



47 Tuc today 
Orbit 

|x a cos(8) = 7.05 mas/yr 
King radius 




-10 



-20 



-30 -40 
I [deg] 



-50 



-60 



-70 



Figure 5. Same as for the top panel of Figure [4] except 
with /i a cos6 = 4.23 mas yr" 1 (top panel) and /i a cos8 = 
7.05 mas yr^ 1 (bottom panel). Corresponds to Figure 6 of the 
original paper. 
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Figure 6. Same as for the top panel of Figure[4] except with fig = 
— 1.54masyr~ 1 (top panel) and fis = —2.56 mas yr —1 (bottom 
panel). Corresponds to Figure 7 of the original paper. 
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Figure 7. Same as for the top panel of Figure [4] except with 
M = 0.9 X 10 6 M (top panel) and M = 1.1 X 10 6 M (bottom 
panel). Corresponds to Figure 8 of the original paper. 
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ABSTRACT 

The Galactic globular cluster 47 Tucanae (47 Tuc) shows a rare increase in its velocity 
dispersion profile at large radii, indicative of energetic, yet bound, stars at large radii 
dominating the velocity dispersion and, potentially, of ongoing evaporation. Escaping 
stars will form tidal tails, as seen with several Galactic globular clusters, however, 
the tidal tails of 47 Tuc are yet to be uncovered. We model these tails of 47 Tuc 
using the most accurate input data available, with the specific aim of determining 
their locations, as well as the densities of the epicyclic overdensities within the tails. 
The overdensities from our models show an increase of 3-4% above the Galactic back- 
ground and, therefore, should be easily detectable using matched filtering techniques. 
We find that the most influential parameter with regard to both the locations and den- 
sities of the epicyclic overdensities is the Heliocentric distance to the cluster. Hence, 
uncovering these tidal features observationally will contribute greatly to the ongoing 
problem of determining the distance to 47 Tuc, tightly constraining the distance of the 
cluster independent of other methods. Using our streakline method for determining 
the locations of the tidal tails and their overdensities, we show how, in principle, the 
shape and extent of the tidal tails of any Galactic globular cluster can be determined 
without resorting to computationally expensive TV-body simulations. 

Key words: Galaxy: globular clusters: general - Galaxy: globular clusters: individual: 
47 Tucanae - Galaxy: kinematics and dynamics - methods: numerical 



1 INTRODUCTION 

47 Tucanae (47 Tuc) is among the most interesting and best 
studied globular clusters of the Milky Way (MW). It is one 
of the most massive, and exhibits some peculiar properties, 
including a recently detec ted rise in the ve locity disper- 
sion at large cluster radii (|Lane et alj|2010ah . Such a rise 
in the velocity dispersion pr ofile of a globular clu ster was 
first noted some time ago bv lDrukier et all l|l998h in M15. 
While a rise of this nature is intuitively unexpected for a 
velocity dispersion profile in Newtonian dynamics, it arises 
naturally in modified gr avitational theories like MOND (e.g. 
IScarpa fc Falomdl2010l 1. 

In the case of 47 Tuc lLane et al | (|2010bh posed, as an 
alternative explanation for this phenomenon, the existence 
of two separate kinematic populations, a relic of its forma- 
tion. A third explana tion was suggested more recently by 
iKiipper et al.l |2010bh . who found a similar rise in veloc- 
ity dispersion in their Wbody models of star clusters orbit- 
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ing in galactic tidal fields. Their alternate explanation arises 
from the contribution of energetically unbound stars which 
haven't yet escaped the cluster p otential. These "potential 
escapers" of lKiipper et al.l (|2010bl ) are preferentially located 
at larger radii, where they dominate the velocity dispersion. 
Moreover, they are an indicator of ongoing, relaxation driven 
mass loss (evaporation), since a star has to become a po- 
tential escaper before it eventually esc apes from the cluster 
to become a member of th e tidal tails (|Fukushige fc Heggid 
|2000| ; IKiipper et alj|2010bl ) . Dynamically cold tidal tails em- 
anating from 47 Tuc are, therefore, to be expected in this 
scenario. Hence, locating these tails will help to clarify the 
mechanism responsible for rising velocity dispersion profiles 
in globular clusters. 

Tidal streams formed in this way have been observed 
for ot her MW globular cluste r s like Palomar 5 and NGC 
5466 (|Odenkirchen et al. I l2003l : iGrillmair fc Johnson! [2006). 
but no tidal tails associated with 47 Tuc have yet been un- 
covered. However. lLane et al. I 2010a|) foun d 25 stars outside 
the tidal radius given in the lHarrisI (|1996T ) catalogue, possi- 
ble evidence for evaporation. Furthermore, these apparently 
extra-tidal stars precisely follow the trend of increased veloc- 
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ity dispersion as found by IKiipper et al.1 (l2010bl) in A^-body 
computations (see Fig. 11 by Lane et al1l2010a ). 

The lack of prominent extra-tidal features could be 
due to the n o n-line ar m otion of stars w ithin tidal tails. 
iKiipper et ail (|2008l ) and lJust et all (120091 ) found that the 
epicyclic motion of evaporating stars results in periodic over- 
and underdensities within the tidal tails, even if they are 
formed by a constant stream of escaping stars. The period- 
icity is spatial, meaning that the overdensities are regularly 
spaced, at least in the case of a circular galactic orbit, and 
are found in both the leading and trailing tails. Their spatial 
separation equates to the distance through which an escaper 
drifts (relative to the cluster) in one epicyclic period, and de- 
creases slowly as the cluster loses mass. Close to the tidal ra- 
dius the escaping stars are accelerated away from the cluster, 
therefore, this motion leads to a pronounced underdensity 
of stars close to the cluster. Hence, most clusters should not 
exhibit extended tidal features, but rather a steeply decreas- 
ing surface brightness close to the cluster (e.g. see Fig. [IJ. 
The epicyclic overdensities, on the other hand, can reach sur- 
face densities several times higher than the average density 
within the tidal tails, but are typically located many tidal 
radii from the cluster and thus the association with the clus- 
ter may not be obvious. The positions of these overdensities 
can be accurately predicted, however, if the dynamical state 
of the cluste r (i.e. its orbit , locati o n fc mass) is sufficiently 
const rained ()Kupper et al 1 l2010al : IKiipper. Lane fc Heggiel 
I2012I ). 

Since the orbit of 47 Tuc has been determined in a num- 
ber of observational studies, and it has been subject to ex- 
tensive numerical modelling, its dynamical state is thought 
to be well understood. Thus, our knowledge of 47 Tuc pro- 
vides an opportunity for predicting and finding the over- 
densities in its tidal tails and, therefore, constraining the 
dynamical state of 47 Tuc independently from other meth- 
ods. Moreover, the detection of the overdensities in the tails 
of 47 Tuc will shed light on the consistency of our theoretical 
understanding of Newtonian gravity. 

This work is organised as follows: based on the most 
accurate data available (Sec. [2j|, we calculate the orbital 
parameters of 47 Tuc and predict the shape of its tidal 
tails and, most importantly, the positions and extents of 
the epicyclic overdensities (Sec. [3). Finally we give a short 
summary and conclusions in Sec. [4] 



2 INPUT DATA 

For this work we require accurate input data to ensure the 
best quality orbital model (see Table [T] for a summary). 
For this reas on we have chosen t o take our positional infor- 
mation from iMcLaughlin et a l. (2006) who claim a nearly 
complete, uniform sample of all stars within l'.5 (~ 5 core 
radii) of the cluster centre, which the authors calculate to 
be at (RA,dec) = (00 h 24 m 05'!67±a07,-72 o 04'52'.'62±0 / . / 26). 
We also take our Helio centric distance information from 
McLaughl in et al.l ([2006) , who calculate a kinematic dis- 
tance of 4.02 ±0.35 kpc by comparing the velocity dispersion 
of stars within the inner ~ 105" with the proper motions in 
the same region to estimate a distance to the cluster. 

For our prop e r mo tion we use the values derived by 
lAnderson fc King! l|2004h . who used Hubble Space Telescope 
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Figure 1. Surface density profile of 47 Tuc ba sed on data from 
lLane et alj l|2010ah , overlaid with a King profile llKinsJ 19621 . solid 
black curve) and a KKBH profile dKiipper et alj 12010a . solid 
green curve). The radius where the surface density drops to zero 
is 56 pc in the case of the King profile. The cut-off radius of the 
core of the KKBH profile is slightly larger (65 pc). However, the 
power-law component of the KKBH extends out to the last data 
point with a slope of —4.8. 



(HST) data to find the proper motion of 47 Tuc to be 5.64 ± 
0.20 mas yr" 1 in RA and -2.05±0.20 mas yr _1 in dec. These 
data are considered the best available data for the current 
work because the Small Magellanic Cloud is in the same 
line of sight as 47 Tuc, allowing the authors to obtain very 
accurate differential proper motions. Indeed they claim an 
order of magnitude greater accuracy than previous studies. 

We have chosen t o use the ra d ial vel ocity of —16.85 ± 
0.16 kms -1 derived by lLane et al.l (|2010al ) from the largest 
known spectroscopic GC study to date. Their survey of 
~ 2200 individual spectra of stellar members provides the 
most accurately determined radial veloc i ty to d ate. The sur- 
face density profile from the lLane et al.1 l|2010al ) dataset can 
be seen in Fig. [T] Note that the King model fit is excellent 
out to large radius (where the influence from the potential 
escapers becomes important, and apparent) and the King 
tidal radius from the fit is in excellent agreement to that 
give n by th e Harris ( 1996) catalogue. Fitting a KKBH pro- 
file (jKiipper et al.ll2010bl ) to the same data yields a slightly 
larger cut-off radius for the core profile (65 pc) and a power- 
law slope of —4.8 at large cluster radii. As we show later, 
the stars at radii larger than the cut-off radius are not extra- 
tidal. The theoretical tidal radius of 47 Tuc is larger than the 
range of the data (see Sec.O, so the fitted cut-off radii are 
both well inside the tidal sphere. Hence, the surface den- 
sity profile does not show a clear cut-off at about 60 pc, 
but continues with a constant power-law slope. Extra-tidal 
signatures like tidal distortions are only expected at much 
larger radii of a few hundred parsecs (e.g. see Fig. [4]| . 
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To estimate surface densities in the modelled tidal tails 
and their epicyclic overdensities, it is important to have an 
independently derived mass-loss rate for the cluster. Using 
detaile d Monte Carlo models o f the dynamical evolution of 
47 Tuc. lGiersz fc Heggid <|20111 ) found an escape rate of 13.9 
stars per Myr, with an average stellar mass of 0.404 Mq, 
giving a mass loss rate of 5.61 Mq Myr -1 . We have adopted 
their value here. 

Many studies exist with measurements of th e total mass 
of 47 T uc, however, the most rece nt of these, lLane et al.l 
(|2010aT ) and lGiersz fc Heggid (|201ll ). based on observations 
and numerical modelling, respectively, find masses of 1.1 x 
10 6 M and 0.9 x 10 6 M which are in remarkable agreement. 
We have, therefore, chosen 1.0 x 10 6 Mq as the baseline mass 
for our models. 

The Solar gala ctocentric distance we take from 
iGillessen et"afl (|2009h who find R = 8.33 ± 0.35 kpc by 
fitting orbits to stars associated with the Galactic centre. 
Since the distance to the the Galactic centre strongly cor- 
rel ates with the mass o f the central black hole (see Fig. 15 
by IGillessen et al"1l2009l ) these orbital fits tightly constrain 
the distance to the Galactic centre, independent of other 
methodologies. 

In a review of twenty independent measurements 
of the Solar circu l ar ve locity dating back to 1974, 
iKerr fc Lvnden-Belll (|l986l ) concluded that a case can 
be made for accepting the IAU standard value of 
220 km s -1 , since the mean value of these twenty stud- 
ies is 222.2 km s _1 . We have, therefore, chosen to use the 
IAU standard of 220kms~ for the current study (see also 
IKoposov. Rix fc Hogdlioiol . who find V c = 224 ± 13 km s -1 , 
in agreement with the IAU standard). 

For the Local Standard of Rest we are adopting 
(U,V,W)^n = l±U+o5l 12.24±° ;£, 7.25±g;|S) kms' 1 by 
ISchonrich et al.l (|2010h as this is the only study to take into 
account the metallicity gradient in the Galactic disc. With- 
out adjusting for this gradient, the V component is signifi- 
cantly affected due to the correlation between stellar group 
colours and the radial gradients of their properties. 

The coordinates of the Galactic North Pole, which are 
required for c oordinate and prop e r mot ion transformations, 
we take from lReid fc Brunthalerj l|2004h . All input data are 
summarised in Table [1] 



3 RESULTS 

Using the input data described in the previous section, we 
first derive the orbit of 47 Tuc about the Galactic centre by 
means of numerical integrati on. We then use this orbit to ap - 
ply the method described bv lKiipper. Lane fc Heggiel (|2012T l 
for predicting the shape of the tidal tails of the cluster. Fi- 
nally we discuss the epicyclic overdensities within those tails. 



3.1 The orbit 

The last 500 Myr of the (x,y) and (x,z) orbital paths of 47 
Tuc derived from the proper motion and radial velocity of 
the cluster (see Table[T| are shown in Fig. [2] For this purpose 
we integrated the cl uster orbit, assuming the Milky Way po- 
tential suggested by I Allen fc Santillanl (l99ll) . consisting of 
a central point mass, a disc and an isothermal halo. This 




Figure 2. The orbital path of 47 Tuc from 500 Myr ago to the 
present time in the (x,y) plane (top panel) and in the (x,z) plane 
(bottom panel), based on parameters from the literature (see Ta- 
bic [TJ . The present day position of 47 Tuc is marked as a plus 
(+) symbol and the position of the Sun is marked as a cross (x). 



potential has a circular velocity of 220kms _1 at the galac- 
tocentric distance of the Sun. In total, we integrated the 
orbit for 1000 Myr back in time, which is the time needed 
to build up the tidal tails out to large angular distances from 
the cluster centre. 

The derived orbit yields a perigalactic distance, R per i = 
5.90 kpc and an apogalactic distance, R apo = 9.18 kpc. 
Hence, the eccentricity, e, of the orbit is 



Ra 



Rr> 



Rapo ~t~ Rpe 



= 0.22. 



(1) 



The current galactocentric distance, Roc, of 47 Tuc is 
7.59 kpc, and the cluster is on its way to perigalacticon. The 
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Table X. Input parameters 



Pctrciirictcr 


Value 


R,6i6r6nc6 




A 09 Inr 

'i.U^ JVJJL. 




Q 47Tuc {o^,V\J\.J) 


nn h 94 m n c i s (S7 

UU ^1 UJ . u / 


MrTantrhliTi pt al (9nnfi1 


°47Tuc W ZUUU J 


— 72°04 /c i2 // fi2 


MrT.aiiP-hlin pt al (900fi1 

1V±L, J_J<X Ll^llllll L^U til . ^ iUUU / 




5.64 mas yr — 


Anderson & King (2004) 


Ms 


-2.05 masyr - 1 


Anderson fc King (2004) 


v r 


-16.85 kms" 1 


Lane ct al. (2010a) 


rt 


56 pc 


Harris (1996); this work 


Aftot 


1.0 x 1O 6 M 


Lane et al. (2010a); Giersz & Heggie (2011) 


A> 


13.9 Myr" 1 


Giersz fc Heggie (2011) 


M 


5.61 M Myr" 1 


Giersz fc Heggie (2011) 


M/N 


0.404 Mq 


Giersz & Heggie (2011) 


(X, Y, Z) G 


(-8.3,0.0,0.0) kpc 


Gillcssen ct al. (2009) 


V® 


220.0 kms -1 


Kerr & Lvnden-Bell (1986) 


(U, V, W) LS r 


(11.1,12.24,7.25) kms" 1 


Schonrich et al. (2010) 



orbital phase, p or b, as defined bv lKiipper et al] (|2011al ) is, 
therefore, 



Porb — 



Rgc Rgc — Rp 
\Rgc\ Rapo — Rp 



= -0.52, 



(2) 



which means that the cluster is half-way to perigalacti- 
con. Its current orbital velocity is 206 kms" 1 , while it is 
IGOkms" 1 at apogalacticon and 264kms~ 1 at perigalac- 
ticon. Hence, the cluster and its tails are currently being 
strongly accelerated. 



3.2 The tidal tails 

For predicting the shape and extent of the tidal 
tails of 47 Tuc, we apply the method introduced by 
iKiipper, Lane fc Heggie] \20li ). Starting at 1000 Myr in the 
past, we follow the orbit of 47 Tuc and calculate its theoret- 
ical tidal radius, XL(t), at time t for each timestep using 



x L (t) = 



GM 



Q 2 -d 2 $/dR 2 



1/3 



(3) 



where G is the gravitational constant, M is the assumed 
cluster mass, $ is the Galactic potential, f2 is the angu- 
lar velocity of the cluster about the Galactic centre and 
R is its galactocent ric distance (e.g. iHeggie fc Hutl 120031 ; 
IKiipper et aill2010al ). 

Into the combined potential of the Galaxy and the clus- 
ter we inject a stream of test particles, each of which may 
represent an escaping star. At the time of its injection each 
test particle is placed at one of the two Lagrange points, LI 
or L2, which lie at a distance of ±xl from the cluster centre 
along the galactocentric radius vector of the cluster (the way 
in which the initial velocity is chosen is described below). 
Since the gravitational attraction of the cluster equals the 
repulsive forces of the effective Galactic potential at these 
points, the test particles can easily escape into the tidal tails 
of the cluster. We then integrate each released test particle 
to the present time, ignoring the mutual gravitational inter- 
action between the test particles but taking the gravitational 
attraction of the cluster into account. 

Since the tidal radius changes rapidly when the clus- 
ter goes through perigalacticon, many test particles can get 
recaptured on the way to apogalacticon as the tidal radius 



quickly grows to its maximal value at apogalacticon. There- 
fore, it is no longer practical to release our test particles 
from the Lagrange points. Instead we introduce a minimum 
radius from which we release test particles, x e d ge , which we 
denote the edge radius. That is, the test particles are re- 
leased from a radial offset, Ax, from the cluster given by 



Ax(t) = max (xL{t), x e d ge ) ■ 



(4) 



This edge radius is found experimentally by increasing 
this lower limit from the perigalactic tidal radius up to 
the value at which no test p article is recaptured (see also 
IKiipper, Lane fc Heggiell2012l . for a more detailed descrip- 
tion). However, we will continue to use the phrase "Lagrange 
points" when referring to the points from which we release 
test particles. 

From the results of the numerical integrations we obtain 
the positions of the test particles at the present time, giv- 
ing streaklines such as those plotted in Fig. [3] Streaklines 
are a concept from fluid dynamics, in which test particles 
are released into a fluid from a given point to visualise flow. 
In engineering streaklines are often produced with smoke 
or dye for trac king flows of air or liquids, re spectively (see 
Section 2.3 bv lKiipper, Lane~ fc Heggie |2012| . for a detailed 
description of the use of streaklines in the current context) . 
Specific to the current study, we let x(t;£o) be the galacto- 
centric position at time i of a particle which was injected 
at time to- Then a streakline in Fig. [3] is the set of points 
x(t;to) in which t is the present time and h ranges over the 
duration of the simulation. Thus streaklines are not orbits. 
By contrast the red curve in Fig.[3]is the orbit of the cluster, 
the locus of its galactocentric position x(i), in which t ranges 
over the duration of the simulation. Such plots are commo n 
in work on tidal tails, e.g. Fig. 20 bv lEvre fc Binnevl (|201lh . 

Our simulations are designed to locate the overdensities 
in the tidal tails of 47 Tuc caused by the epicyclic motions 
of the escaping stars (see Sec. 13. 3[) . These overdensities are 
the direct result of the motions of the escaping stars. Es- 
caped stars have the lowest relative motion, with respect to 
their neighbours, within the epicyclic loops, leading to stars 
'bunching up' at these locations, resulting in overdensities. 
We performed several simulations to cover uncertainties in 
the cluster mass, proper motion and Heliocentric distance of 
the cluster. Therefore, we will first discuss a reference model, 
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Figure 3. The final timestep (present time) of the reference 
model is shown. The black points are test particles. Due to the 
manner in which they were released from the cluster, these par- 
ticles can be considered streaklines, which allow for a simple vi- 
sualisation of the positions, at the present time, of stars which 
have escaped from 47 Tuc in the recent past. These have been 
produced by releasing test particles from the tidal radius of the 
cluster (see text). The solid red curve is the orbital path of the 
cluster. The current position of the cluster is marked by a green 
cross, and its King radius as derived from its surface density pro- 
file (see Fig. [TJ is given by the green dashed curve. Escaped stars 
are slowest within the epicyclic loops, which is where overdensitics 
will be visible. 
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then we will vary single parameters to study the effect of the 
changes. 

Table Q] describes the parameters we have chosen for 
our reference model. This reference cluster has a total final 
mass of 1.0 x 10 6 Mq and loses one star every 0.075 Myr. 
Its proper motion is given by n a — 5.64 mas yr _1 and 
(is = — 2.05 mas yr -1 and it has a Heliocentric distance of 
4.02 kpc. We run our simulations for 1000 Myr, finishing at 
the present time. We have chosen 1000 Myr as this is long 
enough for the first, second and third order overdensities in 
the tidal tails to become fully populated. 

In Fig. |3] the streaklines of the above setup are shown. 
To produce these lines we released the test particles with 
exactly the angular velocity of the cluster from exactly the 
Lagrange points (Eq. The edge radius for this cluster 
mass was found to be 158 pc, whereas the actual tidal ra- 
dius of the cluster varies between lllpc at perigalacticon 
and 171 pc at apogalacticon. The epicyclic motion of the 
test particles in the tails is obvious. The epicyclic loops are 
strongly influenced and disturbed by the cluster mass. 

In Fig. 0]the simulated tails of 47 Tuc are shown for two 
different sets of escape conditions. For the 'warm' escape 
conditions (top panel), the escaping stars retain the angular 
velocity of the cluster plus they obtain a random offset in 
velocity drawn from a Gaussian distribution with a FWHM 
of lkms -1 . Moreover, a Gaussian offset with a FWHM of 
25% of the tidal radius has been added about the Lagrange 
points so that they do not all escape from a single point 
on the edge of the cluster. For the 'hot' escape conditions 
(bottom panel in Fig. 3} we doubled the above fluctuations, 
i.e. a Gaussian velocity offset with a FWHM of 2kms~ 
and a Gaussian spatial offset with FWHM of 50% of the 
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Figure 4. Same as Fig. [3] except with 'warm' (top panel) and 
'hot' (bottom panel) escape conditions of the test particles (black 
points) as described in the text. Each test particle has been given 
a random spatial offset and a velocity offset at the moment of 
escape. Because of these offsets the test particles no longer lie 
along a smooth curve as in Fig. [3] The larger the spread in es- 
cape conditions, the lower the peak density within the epicyclic 
overdensities, however, the locations of the overdensities are only 
negligibly affected. 



tidal radius. However, from lKupper, Lane fc Heggiei (|2012h 
it is obvious that the scatter in escape conditions is in fact 
not very large. Hence, the 'hot' escape conditions can be 
regarded as the 'worst case scenario'. 

We see that adding random fluctuations to the escape 
conditions of the test particles increases the width of the 
tails and scatters their orbits about the ideal orbits shown 
in Fig. [3] The larger the scatter in escape conditions, the 
smaller the peak density in the epicyclic overdensities. How- 
ever, the positions of the overdensities is not altered signifi- 
cantly. 

Furthermore, we have produced several simulations 
to test how the uncertainty in the input parameter val- 
ues affects the output of our model. For proper motion, 
we produce models with fj, a = 4.23 mas yr _1 and /i a = 
7.05masyr _1 as well as with fig = — 1.54 mas yr -1 and 
fls — —2.56 mas yr _1 . For our decreased and increased mass 
models we employ a cluster mass of 0.9 x 10 6 Mq and 
1.1 x 10 6 Mq, respectively. We also include Heliocentric dis- 
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tances of 4. 7kpc and 3.3 kpc beca use, during the writing of 
this paper, IWoodlev et al.l (|2012h published a Heliocentric 
distance of 4.7 kpc, a discrepancy of about 15% from the 
value we assume for our reference model. For these vari- 
ations we assumed the same setup as for the 'warm' case 
mentioned above, except for the parameter which was ex- 
plicitly altered. See Table [2] for a summary of the positions 
and peak densities of the first and second order epicyclic 
overdensities in the trailing tail resulting from each model. 
These models will be further discussed in the following sec- 
tion. 



3.3 Epicyclic Overdensities 

Despite the density of the epicyclic overdensities changing 
slightly from the warm to the hot escape conditions (see 
Models 1 and 2 in Tabled, their positions change very little, 
therefore, for the remainder of this paper we will mostly 
focus on the warm tidal tails, in the interest of brevity. 

Density contours of the warm tidal tails can be seen 
in Figures to The epicyclic overdensities exhibit sur- 
face densities ~ 3 — 4% above the average Galactic back- 
ground in the region near the tidal tails, depending on the 
model (see Table [2J). Tidal features exhibiting very low sur- 
face brightness variations are r ecoverable using matched fil - 
tering techniques. For example, Grill mair fc Johnson! (|2006t ) 
resolved tidal tails with extremely low surface brightness 
over a region of sky with a mean stellar background of 
~ 4500 deg~ 2 , about 50% higher than the background near 
the tidal overdensities of 47 Tuc. Furthermore, after apply- 
ing an "optimal contrast filte r" (effectively a matched fil- 
ter), [Odenkirchen~^tjjI] l|2003h detected features with a peak 
~ 20% above the background counts, and much smaller de- 
viations from the background between the density peaks, for 
the tidal tails of the globular cluster Pal 5. Indeed, the au- 
thors claim that the density peaks are about twice the mean 
density of the tidal tails overall (i.e. the mean observed den- 
sity of the tails is 10% above the background counts). As- 
suming their claimed 4.3 times increase in contrast using the 
matched filter, the density peaks for the epicyclic overden- 
sities in the tidal tails of 47 Tuc will be > 13% above the 
background, well above the detection threshold. 

Furthermore, it can be seen that altering various param- 
eters also affects the positions of the first and second order 
overdensities differently. This is due to the different influence 
of the cluster mass and the orbital phase on t he appearance 
of the tidal tails and the overdensit ies (see iKiipper et al.l 
l2010al ; IKiipper. Lane fc Heggiel I2012T I. Hence, not only the 
location of the overdensities with respect to the cluster is of 
interest when constraining the dynamical state of the clus- 
ter, but also the relative distance between the overdensities 
themselves. 



3.4 Comparison with TV-body simulations 

Figures 7-10 by IKiipper. Lane fc Heggiel (|2012h show that 
our formalism for releasing test particles into the tidal tails 
yields very good agreements with results from direct TV- 
body simulations. However, in comparison to the models 
presented in that study, 47 Tuc moves on a slightly more 
complex orbit which also includes periodic disc crossings. 
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Figure 5. Contour maps of the tidal tails of 47 Tuc as shown in 
Fig. [4] In the top panel the warm escape conditions are shown, 
whereas in the bottom panel the test particles were released with 
hot escape conditions. The contour lines refer to 15, 30, 45, 
stars/deg 2 . The red solid curve is the orbital path of the cluster. 
The current position of the cluster is marked by a green cross, 
and its King radius as derived from its surface density profile (see 
Fig. [TJ is given by the green dashed curve. As can be seen in 
Fig. [4] the peak density in the epicyclic overdensities gets much 
lower when the spread in escape conditions increases. 



Depending on where a cluster crosses the Galactic disc, disc 
shocks can significantly he at up the cluster and incre ase its 
overall mass-loss rate (e.g. IVesperini fc Heggiel Il997). Fur- 
thermore, Eq. [3] which we use for calculating the distance 
from which we release the test particles, has to be taken as 
a rough approximation during these events. Therefore, we 
have to check if our assumptions are still valid under these 
conditions. 

Unfortunately, a direct TV-body simulation of 47 Tuc is 
currently not feasible due to the large number of member 
stars (about 2 x 10 6 ). Instead, we set up a moderately high 
resolution toy model of 47 Tuc with 1.1 x 10 s particles us- 
ing th e publicly available code McLusteeQ (Ki ipper et al.l 
l2011bl ). We used a spectrum of particle masses, m, rang- 
ing from 4 — 40 Mq and a mass distribution proportional to 



www . astro . uni-bonn . de/~akuepper/mcluster/mcluster . html 
www. astro .uni-bonn. de/~webaiub/german/downloads .php 
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Table 2. The density and Galactic coordinates of the first and second order overdensities, the proper motion, final (present day) cluster 
mass, Heliocentric distance and escape conditions for each of our models. Only the trailing tail is represented here for brevity. Information 
is for the peak density in each overdensity only. Note that the approximate average density of the stellar background in the regions of 
the first and second tidal overdensities is 3000 stars per square degree, based on 2MASS point sources. 
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Figure 6. Same as for the top panel of Fig. \E\ except with 
fi a = 4.23masyr _1 (top panel) and fi a = 7.05masyr (bot- 
tom panel). 
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Figure 7. Same as for the top panel of Fig. \E\ except with fig = 
— 1.54masyr — 1 (top panel) and fig = — 2.56masyr~ 1 (bottom 
panel). 



m , adding up to a tota l mass of l.OxlO 6 M©. The par- 
ticles were distributed in a IPlummerl l|l91ll ) sphere with a 
half-mass radius of 10 pc, largely consistent with the surface 
density profile of 47 Tuc. The cluster was in tegrated along 
the orbit of 47 Tuc for 1 Gyr using Nbody6 (|Aarsetbll2003T ) 
on a GPU computer at AlfA Bonn. Over the course of the 
simulation the cluster model lost about 3% of its mass into 



the tidal tails such that its final mass resembles the present- 
day mass of 47 Tuc, within observational uncertainties. 

In Fig. [10] the cluster and its tails are shown after the 
integration. Given the simplicity of our Af-body model, the 
distribution of particles in the tails resembles our predictions 
from the streakline models very well (c/. Figs. [4] & [5}. The 
overall shape and the positions of the overdensities matches 
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Figure 8. Same as for the top panel of Fig. [5] except with M = 
0.9 X 10 6 Mq (top panel) and M = 1.1 X 10 6 M (bottom panel). 
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Figure 9. Same as for the top panel of Fig. [5] except with 
Dhelio = 3.3kpc (top panel) and DHeiio = 4.7kpc (bottom 
panel). 



the predictions from the 'hot' escape conditions best. How- 
ever, a more realistic simulation with ~ 20 times the num- 
ber of particles would have a longer relaxation time and, 
therefore, its mass-loss rate would differ from that of our 
model. Hence, the particle density in this simulation can- 
not be taken as direct proxy for the expected stellar density. 
Moreover, a simulation with 2 x 10 6 particles would exhibit 
much less scatter in escape co nditions and, theref ore, dy- 
namically colder tidal tails (see iKiipper et al]|2008l ). Thus, 
we expect that the tidal tails of 47 Tuc rather resemble the 
'warm' escape conditions. 

Our simplified model (Section I3.2|l assumes that the 
injection rate is constant, however, in the iV-body model 
described in the current Section, the Galactic potential in- 
cludes a disc component. This may, in principle, modulate 
the escape rate via disc shocking. However, this simulation 
demonstrates that disc shocks do not significantly affect the 
shape of the tidal tails or the positions of the epicyclic over- 
densities within them. Furthermore, while it is true that 
the additional energy input from tidal shocks, which is dis- 
tributed within the cluster, r esults in a higher ave rage mass 
loss rate, it was estimated bv lAguilar et alj |l988T ) that such 
shocks only contribute ~ 1% to the total destruction rate. 
Moreover, this increase in mass loss is distributed nearly ho- 
mogeneously over the whole orbit because the escape of the 
stars is delayed until well after the disc shock, often for sev- 



eral orbits of the st ar within the cluster, even for clusters on 
non-circular orbits (jKiipper et al.ll2010al ). This delay is the 
result of the stars preferentially leaving the system through 
the Lagrange points, where the potential barrier is at its 
minimum. Stars excited to escape energies by shocks leave 
the cluster on about the same time-scale as due to other 
evaporative processes, and may even stay bound to the clus - 
ter for as long as a Hubble time (|Fukushige fe Heggid[200ol ) , 
at least for clusters on circular orbits. The energetically un- 
boun d stars enhancing the outer velocity dispersion of 47 
Tuc (|Lane et al.ll2010aP ) also support the concept of delayed 
escape, especially if these unbound stars were excited during 
the last disc passage which, from our simulations, occurred 
about 70 Myr ago. Therefore, because energetically unbound 
stars cannot all escape immediately following shocks, and 
can stay bound to the cluster for long periods, we do not ex- 
pect to see additional substructure in the tidal tails caused 
by disc crossings. 

While we us e d the analytic potential suggested by 
lAllen fc Santillanl (|l99ll ) f or the streakline models as well 
as for the ./V-body model, iDehnen fc Binnevl (| 19981 ) argue 
that the Galactic disc structure is actually more complex 
than the single-component exponential disc used by Allen 
& Santillan. However, the strongest disc shock 47 Tuc ex- 
periences is at 6 kpc. At this distance the two potentials do 
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Figure 10. Distribution of tail stars from a simple TV-body simu- 
lation. The shape of the tails and the positions of the overdensitics 
agree very well with the predictions from our streakline models. 
A more realistic TV-body model with a larger number of particles 
(here we use 110 000) will show less scatter in escape conditions 
and will result in dynamically colder tidal tails (see Sec. 13.41 1. Par- 
ticles within the King radius are not shown in the upper panel 
for clarity. The contour lines in the lower panel refer to 2.5, 5, 
7.5, particles per square degree. Due to the smaller number of 
particles in the tidal tails we used a Gaussian smoothing kernel 
with FWHM = 1°, twice as large as for our other contour plots 

(Figs. EE}. 

not differ greatly (e.g. iKaempf. de Boer fc Altmannl 120051 : 
IJilkova et afl|2012h . Therefore, we do not expect our predic- 
tions to be signif icantly affected by t he str onger disc shock- 
ing suggested by iDehnen fc Binnevl ll 19981 ) . 



4 CONCLUSIONS 

We have employed the most accurate available data on the 
Milky Way globular cluster 47 Tuc to compute its orbit 
about the Galaxy. Moreover, w e used the streakline method 
l|Kiipper. Lane fc Hcggic 2012) to predict the shape of the 
tidal tails of 47 Tuc and the locations of the epicyclic over- 
densities within those tails. By varying the input parameters 
of our models we show how the most important input val- 
ues influence the locations and peak surface densities of the 
epicyclic overdensities. 



We find that the epicyclic overdensities exhibit peak 
surface densities ~ 3 — 4% above the Galactic background. 
These should be observationally detectable, assuming a 
matched filter, or similar analysis technique, is used. 

Clearly variations in Heliocentric distance and proper 
motion alter the positions of the tidal overdensities more 
than variations in any other parameter in our model. Our 
variations in proper moti on are much larger than the ob- 
servational error stated by I Anderson fc Kind |2004l ) and we 
find only a small change in the position of the overdensities. 
However, determining the distance to 47 Tuc is an ongoing 
problem, with many different e stimates published over many 
years (see lGratton et aL 1 12001 for several examples). There- 
fore, finding these overdensities observationally will not only 
provide further constraints on the proper motion of the clus- 
ter, but will also constrain the other dynamical characteris- 
tics, including, importantly, the Heliocentric distance. 

fn addition, uncovering these overdensities observation- 
ally may reinforce our understanding of Newtonian gravity. 
While it is not within the scope of the current paper to make 
predictions regarding the tidal tails and their overdensities 
in modified gravity environments, it would be extremely in- 
teresting to see the morphology of the tidal tails as predicted 
by modified gravitational theories, since we would expect 
them to differ significantly from our predictions. 

Finally, our study has shown how, in principle, the 
shape and extent of the tidal tails of any star cluster can 
be quickly predicted without requiring extensive A-body 
simulations. This will be especially helpful for large param- 
eter studies or for massive star clusters which cannot be 
addressed by means of direct A-body computations due to 
hardware limitations. 
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